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Abstract 

We consider Non-Standard neutrino Interactions (NSI), described by four-fermion operators of the form (^q7^/3)(/7/), 
where / is an electron or first generation quark. We assume these operators are generated at dimension > 8, so 
the related vertices involving charged leptons, obtained by an SU(2) transformation us — ► eg, do not appear at 
tree level. These related vertices necessarily arise at one loop, via W exchange. We catalogue current constraints 
from sin 2 8w measurements in neutrino scattering, from atmospheric neutrino observations, from LEP, and from 
bounds on the related charged lepton operators. We estimate future bounds from comparing KamLAND and solar 
neutrino data, and from measuring sin 2 8w at the near detector of a neutrino factory. Operators constructed with 
Vfj, and Ve. should not confuse the determination of oscillation parameters at a j/factory, because the processes we 
consider are more sensitive than oscillations at the far detector. For operators involving v T , we estimate similar 
sensitivities at the near and far detector. 



1 Introduction 



Lepton flavour violation is observed in atmospheric and solar neutrino experiments. The atmospheric neutrino deficit 
is mainly due to muon neutrino disappearance: the flux measured at Super-Kamiokande (SK) has a strong zenith 
angle dependence, which deviates from the Standard Model (SM) expectation by more than la pQ. Neutrinos change 
flavors as they travel to the Earth from the center of the Sun, as was seen directly at the Sudbury Neutrino Observatory 
(SNO) through the measurements of the charged current and the neutral current reactions for S B solar neutrinos [2]. 

These results demonstrate simply that new physics is required. In principle, the atmospheric and solar neutrino 
deficits can be explained by neutrino masses or by giving the neutrinos new interactions. Data disfavor non standard 
interactions ( NSI) of neutrinos as an explanation for the atmospheric neutrino deficit [3J 31 [5] through the energy and 
baseline dependence. KamLAND detector has recently confirmed the large mixing angle (LMA) oscillation explanation 
of the solar neutrino puzzle [BJ. Prior to this, NSI were a viable alternative solution pj [5] El- 
Neutrino masses are the leading mechanism in the solar and atmospheric anomalies. Nonetheless, NSI may be 
comparable to (or larger than) oscillation effects in other processes or at other energies. This is a particularly relevant 
issue for neutrino factories [TUl HU [12 [13j [HI [15l [TBI [17] , where NSI may affect the oscillation parameters inferred 
from experimental data, biasing the value of some of the mixing angles. For instance, it was suggested in [16] that NSI 



* Reference [9] found that NSI induced in the R-parity violating MSSM could not explain the solar neutrino deficit. However, their 
operators were otherwise constrained to be at least an order of magnitude smaller than the solutions found in [8] . 



could interfere with the determination of sin 2 #13 at a neutrino factory!]]. It is important to understand how significant 
is this possibility. The original aim of this paper was to show that short baseline and precision experiments are more 
sensitive to NSI than the far detector at a neutrino factory. We will see that this is true for some NSI, but borderline 
for those involving a v r . 

We consider neutral current NSI, from a phenomenological perspective: we follow |18j . and assume that the new 
physics which induces the non-standard (pp)(f '/) operator, where / is a charged lepton (£) or quark, does not introduce 
new charged lepton physics at tree level. Within this approach, NSI can be constrained [TH [19] from neutrino deep 
inelastic scattering experiments and from elastic scattering v — e, where the baseline is too short for oscillations. 
NSI would contribute to v scattering events, and therefore to the determination of sin 2 9yy in these experiments. So 
short baseline, high- flux neutrino experiments, that measure for instance sin 2 9w, can set significant bounds on NSI 
involving v e and v^. 

A {vv)(f f) operator, will nonetheless induce a (tt)(j '/) and/or (tv)(f '/') operator, where /' is the SU(2) partner 
of /, via external one-loop Standard Model dressing. This is independent of the new physics that induces the [vv)(f f) 
operator. Even in the case in which there are no tree level bounds on NSI of neutrinos, because there is no appropriate 
experiment, in general radiative corrections will generate other types of interactions which could be tested at present 
(or future) experiments. Bounds on charged lepton operators therefore set model independent bounds on the NSI 
operators. However, at present these bounds are only significant for /1 — e flavor changing operators, because the loop 
suppression is a small number. NSI in loops can also make flavour dependent contributions to the decay rates of the 
electroweak gauge bosons, which can set relevant bounds on flavour diagonal NSI. 

Present bounds still allow sizable (diagonal) NSI of v T , but we will show that they can be constrained by comparing 
SNO/SK solar data and KamLAND results. The bounds are significantly better if KamLAND finds an oscillatory 
signal, that is if Am^, ^ 10 -4 eV, as we will assume for definiteness. 

We introduce our notation and assumptions in section 2. Section 3 presents current constraints on flavour diagonal 
and flavour changing neutral current NSI, both from tree level effects (short baseline, high-flux neutrino scattering 
experiments, LEP, atmospheric neutrinos) and from one-loop processes (flavour changing charged lepton interactions, 
LEP). In section 4, we discuss the sensitivity of solar neutrino experiments, SNO and Super-Kamiokande, using 
KamLAND to pin-point the solar oscillation parameters, Am 2 o/ and 6*12. We also discuss the sensitivity of future 
experiments, using the near detector at a nufactory as an example. We summarize our results in section [SJ tabulating 
the best current and future bounds we obtained. 



2 Notation and Assumptions 

At energy scales <§C raw (where there is a large amount of precise v scattering data) , the Standard Model interactions 
of neutrinos can be described by the effective Lagrangian 

C eff = -2V2G F {[D plp L£p][fyLf'] + h.c.) - 2^2G F £ 9 P [^l P LM[h P Pf] (1) 

P,f,P 

where P = {L, R}, £ is a. charged lepton, / is a lepton or quark, /' its SU(2) partner, and the Z couplings g P are given 
in table [TJ Greek indices from the beginning of the alphabet (a,/3,...) label lepton flavours, roman indices 
correspond to neutrino mass eigenstates, and late alphabet Greek letters (p,a...) are space-time indices. 

We consider non-standard, neutral current neutrino interactions, so we add operators with the form of the second 
term in equation (JTJ) . We do not include new charged current interactions. As discussed in jTOl [TU [TJ] [T5], [2D] , NSI can 
contribute to a " neutrino oscillation" signal via charged current interactions in the source or detector, or via neutral 
current interactions in the propagation from source to detector. However, we anticipate that other experimental 
processes are more sensitive to charged current NSI that long baseline neutrino oscillations (e.g. flavour changing NSI 
in the source could induce taus and wrong sign muons in the near detector). 

Non-renormalisable operators involving a Standard Model neutrino and anti-neutrino can be ordered by their 
dimension, or by their number of legs in the SU(3) x U(l) invariant effective theory of SM fermions and photons. 
These options are different, because Higgs fields saturated by the vacuum expectation value (vev) are not counted as 
legs. If we count by legs, then new physics coupled to neutrinos can appear as a four fermion operator. We require 
this operator to conserve electric charge, colour and lepton number, which forces it to be of V ± A form. We assume 
three light neutrinos, with Majorana masses and the SM interactions of equation (UJ). We allow them to have NSI, 
parametrised as 

C eff = -e{p^/2G F {v alp Lv )Ch P Pf) (2) 

t This interference might be solved through the tail of the spectra at large energies (Ref. I17p . 
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Z couplings 


f f 
91 9r 


e,fj,,r 
u, c, t 
d, s, b 


\ o 

— ^ + sin 2 6 W sin 2 6>w 

1 2*2/) 2*2/) 

2 — 3 sm uw — g sin a\v 
— 5 + 3 sin 2 9w ^ sin 2 9w 



Table 1: Z couplings to SM fermions. 



where / is now a first generation SM fermion: e, u or d, and P = L or i?. We are not concerned with / from the 
second or third generation, because such interactions could not affect oscillation experiments. We neglect possible CP 
violation in the new interactions (this has been considered in [T3l[2T]), so we take s j^ € 5ft. If the only new neutrino 
interactions are neutral current, the neutrino flavour basis is well-defined (see [20), and we can label neutrinos by 
their SM flavour. 

Assuming that the four-fermion vertices of equation ([2]) arise in an SU(2)®U(1) gauge invariant theory containing 
the SM spectrum with a single Higgs doublet, they can be generated by operators of dimension six, eight and larger 
[18j containing more and more vevs of the Higgs doublet. In this paper we will not consider dimension seven AL =2 
operators [22j because, unless their coupling is very small, we expect these interactions to generate too large neutrino 
Majorana masses. There is an important difference between the dimension 6 and dimension 8 (yv)(f f) operators: new 
physics which induces the dimension 6 operator also induces an operator involving charged leptons, with a coefficient 
of the same order (by SU(2) invariance) [23 . Charged lepton physics imposes tight constraints on these coefficients 
of dimension 6 operators. At dimension 8, an operator as in equation (J2J) can appear at tree level, for instance as 

e R (H^a a £){£a a H)e R -> -\{H) 2 'i&f Re)^ > p Lv) (3) 

without any charged lepton counterpart [THj 0. We include the operator and neglect the associated tree-level 
charged current vertex, so we are effectively assuming that these operators are generated at dimension eight. 

If new physics operators are generated at tree level by exchange of particles with mass A, naive power counting 
tells us that dimension six operators should give rise to e ~ h 2 v 2 F /A 2 where h is some generic coupling of the new 
physics particles. Taking into account present data collected at LEP and Tevatron one can reasonably assume that 
A > 200 GeV. Even though somehow lighter particles are not excluded (e.g. if neutral or only produced in pairs or 
having small couplings with ordinary particles), we are using an effective Lagrangian approach which is only valid 
below A. This counting tells us that we expect e < h 2 for dimension 6 operators. If NSI interactions are generated 
by dimension 8 operators one expects an extra suppression proportional to v F /A 2 which is only important if A ^> Vp, 
therefore one naturally expects e < 1. However, h could be relatively large without leaving the perturbative regime. 
Thus, e's order one are not completely unnatural if the scale of new physics is not extremely large and the couplings 
of new particles are large. It is, therefore, important to check how large these NSI can be on purely phenomenological 
grounds. Thus, in this paper, we will concentrate on the phenomenology of U(l) x SU(3) invariant (z^ Q ;/g)(/ /) 
operators, with arbitrary coefficients. As mentioned in the introduction, there will nonetheless be constraints from 
charged lepton physics on these operators. One loop SM dressing of the tree level NSI vertices (v a v$)(j T f) (for instance 
exchanging a W between the external legs), will necessarily induce (e Q e^)(//) and/or (v a ep)(f /'). This is discussed 
in section I3~2l 

In the effective theory described by equations ([T]) and @ , there will normally be a number of four-fermion operators 
which can contribute to a process. Experimental measurements of different processes will set limits on different sums 
of operators. It is common, in setting a limit on a given four-fermion operator, to assume that all other four-fermion 
operators are zero. This approach makes sense when the new physics contributions add incoherently, e'est a dire when 
one adds probabilities not amplitudes. It also makes sense if the different new physics amplitudes interfere with 
different SM amplitudes — for instance, if there is an NSI - SM interference term cx e q ^ a g q L + s^ t g q R , these two terms can 
only cancel against each other if the NSI know sin 2 0w- We consider this unlikely (since we assume they are generated 
at dimension 8, rather than 6), so we quote limits on one e at a time in our summary tables. 

We will assume that the only unknown new physics are the non-standard neutrino interactions. We calculate 
what current experiments should measure according to the SM, and require the NSI contribution to be less than the 
experimental error (1.6 a), or less than the theory - experiment discrepancy. 

tThe Lorentz indices in eq. ^ are contracted leji. The as carry gauge SU(2) indices. 



3 Present bounds 



Before considering particle physics bounds, let us discuss briefly astrophysical and cosmological bounds on neutrino 
NSI. Neutrino interactions with matter, electrons and light quarks, can affect many astrophysical and cosmological 
scenarios. They could keep neutrinos in thermal equilibrium with ordinary matter for a longer time at the time 
of nucleosynthesis and disturb one of the great successes of modern cosmology. They could also produce stronger 
interactions of neutrinos with matter in the core of supernovae therefore keeping neutrinos trapped for a longer time 
and disturbing the duration of the neutrino pulse observed in SN1987A [24] [25l [26]. They could also contribute to 
the energy loss of stars due to processes like plasmon decay (7* — > vv) which are determinant for the evolution of red 
giants. All these processes occur in the SM mediated by neutral or charged currents and the SM value is essential 
to understand the three scenarios mentioned. Now the question is how large NSI are allowed in order not to disturb 
present observations? For instance, the case of plasmon decay in red giants has been used to place stringent bounds 
on a possible neutrino magnetic moment [24] . The SM gives contributions to plasmon decay from neutral and charged 
currents for v e and only from neutral currents for v^ and v T . To destabilize the SM results one should modify some of 
these interactions by a factor larger than one [23]. Therefore, generically one can say that e < 0(1) from energy loss 
in red giants. We expect, at best, similar results from supernova data or nucleosynthesis. As we will see, laboratory 
data already place better limits on v e and NSI, and SNO, Super-Kamiokande and KamLAND will also improve 
the bounds for v r . 

We first consider tree level effects of the operators in (J2J, which contain only the neutrino current with either the 
electron or first generation quark currents. Low energy scattering experiments can constrain NSI involving v e and v^, 
while to derive bounds on diagonal v T NSI one should use the measurement of the e + e~ — * z/z/7 cross section at LEP 
[18] and atmospheric neutrino data [5]. We set further bounds using the fact that such operators always induce one 
loop effects in much better tested charged lepton processes. These constrain flavor changing NSI involving first and 
second generation neutrinos to be undetectably small. 



3.1 Tree level effects 
3.1.1 v scattering experiments 

Neutrino NSI with either electrons or first generation quarks can be constrained by low energy scattering data. We 
review previous analysis [19] 118] and update them by including the recent results of the NuTeV experiment. As we 
shall see below, the bounds are rather stringent for z/ M interactions, looser for v e and do not exist for (diagonal) v r . 

We present bounds assuming that only one operator is present at a time, for the reasons explained above, although 
we also comment on how these limits are relaxed when several diagonal NSI are considered simultaneously. 

• v & e — > ve scattering 

In the presence of neutral current neutrino NSI the v e e elastic cross section is given by 



. , 2G 2 F m e E v 
<r(v e e — > ve) = 



(1 + 9l + e? e ? + £ K L e \ 2 + -Ag% + et«? + \ £ |4*| 2 



3 



(4) 



where g e L = —0.2718 and g e R = 0.2326 are the SM neutral current couplings of the electron, including electroweak 
radiative corrections and corresponding to the best fit point of the latest SM global fit of precision observables 
(without including NuTeV). 

The most accurate measurement of this cross section is the LSND result [27] : 

a{v e e -> ve) = (1.17 ± 0.17) G ^ m ' E " , (5) 

which, taking into account the SM prediction a(v e e — ► v e e)\sM — 1-0967 G F m e E v /ir , translates into the following 
90% CL bounds on diagonal v e e NSI (assuming only one operator at a time): 

- 0.07 < ef e < 0.11 (6) 
-1. < elf < 0.5 (7) 

We can also set bounds on flavour changing NSI. These are only relevant for v T v e interactions, because for v^v e 
better bounds are obtained from the one loop effects discussed in section I3T21 Assuming there are only flavor 
changing NSI we obtain: 

|4el<0.4 \ef e \<0.7 (8) 
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Figure 1: Bounds on flavor conserving non-standard v e e interactions from LSND experiment. Allowed regions at 90% 
CL are between the two ellipses. 



One can wonder how these bounds would be relaxed when allowing for several operators to be present simulta- 
neously. We consider then both, left- and right-handed diagonal NSI. The result is shown in Fig. [1] The 90% 
CL allowed region is between the two ellipses, and corresponds to the range 



0.445 < (0.7282 + e 



eL\2 



1 



(0.2326 + £ee) < °- 725 



(9) 



• v e q — > vq scattering 

The CHARM collaboration measured the following combination of v e N and v e N cross sections 

a{v e N -> vX) + a{v e N -> vX) 



a(v e N -> eX) + a{D e N -> eX) 



(he) 2 + (~9Re) = 0.406 ±0.140 



(10) 



Since charged current NSI are strongly constrained, we neglect them and use this measurement to bound the 
neutrino neutral current NSI. In this case the effective couplings (gLe) 2 , (<?.Re) 2 are given by 



(gt 



.dL\2 



5>1 



(~9Le) 2 = (SL+^ + EKel 2 

(he) 2 = (9l + e u e ?f + J2\e u «\ 2 + (9t + e d *f + J2\e 



dL\2 



dR\2 
ae I 



(ii) 

(12) 



The SM couplings corresponding to the best fit are (he) 2 = 0.3042 and (he) 2 — 0.0301. Using this result, the 
90% CL bounds on flavour diagonal NSI are 



- 1. < e e " e L < 0.3 
-0.3 < ef% < 0.3 
-0.4 < e** < 0.7 
-0.6 < eff < 0.5 

when assuming only one operator at a time. 

The corresponding 90% CL bounds for flavour changing NSI interactions are 

|4f|<0.5 q = u,d P = L,R 



(13) 
(14) 
(15) 
(16) 



(17) 



Again, these bounds are only relevant for v e v Tl since for ^ e f M tighter ones are derived from one loop effects. 

If we consider all kind of diagonal NSI, the allowed regions at 90% CL are limited by two four dimensional 
ellipsoids and are given by 

0.176 < (0.3493 + e u ^ f + (-0.4269 + e^) 2 + (-0.1551 + e^f + (0.0775 + e d *f < 0.636 (18) 



The CHARM II collaboration gives the following results for vector and axial- vector e — couplings [55] : 

g e, = -0.035 ±0.017 and g\ = -0.503 ± 0.017 (19) 

where they have used LEP forward-backward asymmetry to determine the signs. From these one gets 

g% = -0.269 ±0.017 and g% = 0.234 ± 0.017 (20) 

The SM values of the left- and right-handed couplings are the same as for v e e scattering, so we can readily derive 
the 90% CL bounds on diagonal NSI 

- 0.025 < ef^ < 0.03 (21) 
-0.027 < ef. < 0.03 (22) 



as well as on flavour changing operators, 



when we allow only flavor changing NSI. 



Kl<0.1 P = L,R (23) 



where 



• v^q -> vq 

The NuTeV collaboration measures the ratios of neutral current to charged current neutrino-nucleon cross 
sections, which for an isoscalar target and at leading order are given by 

* - ^^-<W + ;<W. <*> 

ajPN -» fiX) 

a(uN -► nX) ■ [ ' 

Neglecting charged current NSI, the effective couplings (gi M ) 2 , (gn^) 2 are as given in eas. (fTT|) . (fT2|) . just changing 
e — > n in the coefficients of the neutral current NSI. The values of these couplings reported by NuTeV are [3"0] 

(g Ltl ) 2 = 0.3005 ± 0.0014 and {g Rfi ) 2 = 0.0310 ± 0.0011 . (27) 

While {gR^) 2 is in agreement with the SM, (g^)^ = 0.0301, (<7l/J 2 is about 3a away from the SM prediction, 
(Sl^sm = 0-3042. 

The NuTeV result for {g Lli f can be fitted (at 90% CL) by 

- 0.009 < < -0.003 or 0.002 < ef£ < 0.008 , (28) 

and, in principle, since the measured (g^) 2 is smaller than the SM prediction, pure left-handed flavour changing 
NSI are excluded because they do not interfere with the SM amplitude and therefore always give a positive 
contribution. 

Alternatively one can assume some other explanation [31] of the discrepancy, and estimate that NSI should 
contribute less than 1.64 a to the NuTeV result. This leads to the constraints: 



|e«£|< 0.003 q = u,d (29) 

-qLl 

For diagonal right-handed NSI the 90% CL allowed ranges are 



|4jjl < 0.05 q = u,d (30) 



0.008 < < 0.003 (31) 



0.008 < e d * < 0.015 , (32) 




Figure 2: Bounds on flavor conserving non-standard v^q interactions from NuTeV experiment. The 90% CL allowed 
regions are the thick ellipse in the left panel and the region between the two ellipses in the right one. 



while for flavour changing interactions, 

< 0.05. q = u,d (33) 

These bounds are relaxed if we consider several operators present simultaneously. In Fig. [2] we plot the 90% CL 
limits on diagonal neutrino NSI with u- and d-type quarks, both for left- and right-handed interactions. The 
allowed regions are the thick ellipse in the left panels and the region between the two ellipses in the right one. 
They are described by the equations: 

0.2982 < (0.3493 + e^f + (-0.4269 + e^) 2 < 0.3028 (34) 
0.0292 < (-0.1551 + £"£) 2 + (0.0775 + e d S) 2 < 0.0328 (35) 



3.1.2 LEP 



The authors of Ref. [T5] have pointed out the importance of the e + e~ — > vv^j cross section measured at LEP II in 
order to constraint neutrino NSI. For the case of diagonal v r e interactions these are the only laboratory bounds, and 
for v e e they are comparable to the LSND limits already discussed. We refer the reader to [TS] for the detailed analysis, 
and just summarize here the 90% CL bounds when only v T e NSI are considered, which can be read from their Fig. 4: 

E 

- 0.6 < e% L T < 0.4 (36) 
-0.4 < e£? < 0.6 (37) 



This reaction is also useful to constrain flavor changing NSI, but the bounds are comparable or looser than the 
ones derived from elastic scattering, namely 

|e^|<0.4 P = L,R, a = r, = e,n (38) 



§There are indeed two ellipses, but they can not be distinguished in the figure. 
^We have estimated a ~ 10% accuracy at 90% CL. 
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assuming only flavor changing non-standard operators. 



This process cannot constrain v T q interactions. However, it is obvious that any NSI interaction can contribute 
to Z — ► ffvv at LEP1 or ee + — > ffvv at LEP2. These processes occur in the SM with a pair of virtual Z, W, 7. 
They have also been observed and the observations roughly agree with the SM expectations. For instance the cross 
section with virtual Z and 7* going to quarks and neutrinos e + e~ — > Z7* — * gqi/i/ has been measured by DELPHI 
at energies above 189 GeV to be oz~i* — (0.129 ± 0.038) which has to be compared with the SM prediction of about 
0.092-0.084. This suggest that NSI should be at most as strong as SM interactions. It is difficult to extract more 
precise information on possible NSI from these data because it is based on a particular pole structure which is not 
shared by the NSI. Clearly, to be more precise on the NSI a dedicated study should be performed. In addition, NSI 
of v T with quarks can in principle be constrained through one loop effects both in the invisible and the hadronic Z 
width, but typically these bounds are 0(1) (see sect. 13.2( 1. 

One can ask whether TESLA would set improved bounds. The sensitivity at TESLA to contact interactions of 
the form rj2\/2GF{e'y p Pe){f^p^,,e) can be estimated [33] to be ij 0.5 x 10 -3 . Using c ~ .002 from section I3T21 
r) ~ ce implies e < O(0.3). This is not particularly significant, and corresponds to a new physics mass scale below 
500 GeV (y^s at TESLA) so the contact interaction analysis is not appropriate. A more promising channel might be 
ee — * j(vp), but to our knowledge, the limits that could be extracted from searching for this at TESLA have not been 
studied. 

3.1.3 Atmospheric neutrinos 

Exotic (i.e., no-oscillation) solutions of the atmospheric neutrino problem are disfavored by the energy and baseline 
dependence [3J 0]. In particular, in Ref.[5] an explanation of atmospheric neutrino data in terms of pure NSI of 
neutrinos is ruled out at 99% CL, but a combined analysis is performed which includes both, oscillations and non- 
standard neutrino-matter interactions, and allows to set significant bounds on diagonal as well as flavour changing 
neutrino's NSI. Two comments are in order. First, note that such bounds apply to the vector coupling constant of the 
NSI, — + £^0, since it is the only one which appears in neutrino propagation through matter. Second, the 
study does not consider NSI in the production neither in the detection of the neutrinos. However both processes take 
place through charged current interactions, which are better constrained, so one does not expect sizable effects. 

Notice that when NSI appear only in the propagation of neutrinos in matter, the relevant parameters are the 
flavour changing NSI couplings and the difference between the diagonal ones, denoted in [5j[33] as the non-universality 
parameter e'. The combined analysis leads to the following 90% CL bounds when both diagonal and flavour changing 
NSI are simultaneously present [33] : 

- 0.016 < ef T < -0.009 < e% < 0.013 (39) 

\e%-e%\<0m (40) 

Although the fit has been done assuming only neutrino-down quark NSI similar results can be expected for the 
neutrino-electron and neutrino-up quark cases. 

This result leads to a stringent constraint on flavour changing four-fermion operators involving v^v T , and provides 
a complementary bound on diagonal NSI of v T with quarks, only loosely constrained by LEP. 

3.2 One loop effects 

On general grounds we expect that interactions in which the v a are replaced by the corresponding leptons will be 
generated by one loop diagrams with virtual W's. Effective interactions, however, are nonrenormalizable and, therefore, 
a precise prescription has to be given in order to estimate these corrections. Our point of view will be that these are 
originated from a more complete theory at scales A 3> mw which is renormalizable (or perhaps finite) and in which 
observables can be computed in terms of a few parameters. 

We illustrate our point with a simple toy model, in which exact loop calculations can be easily done. It is not 
a realistic example of the type of model we wish to constrain (see e.g. [18] for such models), because it induces 
(<37 p Pe) {Vf/YpLvr) and (f^ p Pe) {pelpLvp) simultaneously at dimension 6 with approximately the same coefficient. So 
t appearance at the near detector of a neutrino factory would be more sensitive to this model than neutrino scattering. 

As a guide for the type of calculations we are going to perform, we imagine the standard model extended by a 
singly charged scalar singlet h + , with mass M, which has the following interaction with the standard leptonic doublet, 




Figure 3: One- loop contributions to four-fcrmion interactions in a theory with a charged scalar singlet. 



£, (£ = it 2 I c ) 

C h =f a p(ljp)h- + h.c. , (41) 

where f a p is a coupling antisymmetric in flavour indices. Expanding in the fields 

£h = 2f a0 (VtZelp) h~ + h.c. . (42) 

Exchange of scalars will generate a four-fermion interaction of the type we are considering (for a review of this model 
from the effective Lagrangian point of view see |34|). 

£ = Ifapfjsj^ (eI>L 7 ) (v La e c Lfj ) = ZfafffySJp (epYLes) {D a jpLv y ) . 
We are only interested in NSI with electrons, so we take f^ T = and find that this model gives 

C = -2V2G F ef T (e 7 p Le) {v^ p Lv r ) + h.c. + • • • , (43) 

where 

eL _Uef? e .mw , . 

g 2 M 2 ' 1 ' 

In addition this model also generates, at tree level, interactions of the type {e^ p Lr) (D p j p Lv e ) and (Jl^Lt) (y^pLv^) 
which are not interesting for our discussion and which we have represented by the dots in eq. 



Although at tree level this model does not provide any contribution to t~ — > /i~e~e + it does contain contributions 
at the one-loop level. The only contribution with and v T in virtual states is given by the diagram in Fig. [3] (note 
that in this model there are other diagrams contributing to the process with other types of neutrinos in virtual states). 
In fact the diagram in Fig. [3] can be easily computed and it is finite. In the limit of M\y, M 3> m T it generates the 
following interaction among four charged leptons (after Fierz reordering) 

C = - 2 V2G F s%-^F (^Pj (eYLe) {jl lp Lr) , (45) 
"~ s w \ m w J 

where e e J^ is given by eq. (|44|) and the function F(x) is 

x 

Fix) = ln(x) ~ ln(x) . 

x — 1 

The important thing is that to obtain this logarithmic contribution we do not need to know all the details of the 
complete theory. We can compute it by using the four-fermion effective interaction (Fig. [4j) . 

This diagram yields a logarithmically divergent contribution. However, the divergent part can be unambiguously 
calculated. In the effective field theory language one should add a counter-term (a 4-charged fermion interaction) to 
the Lagrangian in order to absorb this divergence. The coefficient of this new operator will run under the influence 
of the diagram in Fig. 01 Therefore even if at the scale of new physics, A, there is no 4-charged fermion interaction 
at tree level it will be generated through running from the scale A to the electroweak scale. Thus, to compute the 
logarithmic contribution it is enough to compute the divergent part of the diagram in Fig. 0J The finite parts have to 
do with possible non-logarithmic contributions generated at the scale A which cannot be computed without knowing 



Figure 4: One-loop contributions to four-fermion interactions in the effective theory. 



the details of the complete theory and should be absorbed in the initial conditions for the running of the different 
operators. However, if the scale A is much larger than the electroweak scale, one can reasonably assume that the 
logarithm ln(A/mw) dominates completely the result. 

If we apply this point of view to the model we just considered we obtain exactly the same result as in eq. (|45|) with 
F(M 2 /rriyy) replaced by 21n(A/mvr) + K (A), where k(A) takes into account possible non-logarithmic contributions 
at the scale A. So the effective Lagrangian calculation gives the correct answer if we identify A with the mass of 
the charged singlet, M and take k(M) — 0. Note, however, that in the effective field theory calculation we have no 
way to determine k(A), which depends on the details of the complete theory. In fact, in the effective field theory A 
can only loosely be related to the masses of the unknown more complete theory. k(A), somehow, parameterizes all 
these unknown details. The important point is that we expect it to be at most order 1 and negligible in front of the 
calculable logarithmic piece if A is large enough. 

In order to set reliable bounds on the e's using these one loop corrections we would need to know roughly the size 
of A. However this is not known. Bounds will be set directly on the e's at tree level and on e\n(h/mw) at one loop. 
Of course, e also depends implicitly on m^y/A 2 as in eq. (|44[) . however, it also depends on other parameters, in this 
case the Yukawa and gauge couplings, which, at the level of the effective four-fermion theory cannot be completely 
disentangled. One might expect the logarithm to give some enhancement, but, since the size of this enhancement 
cannot be reliably computed we choose to be conservative and take \n(A/mw) w 1 in all the bounds we will set from 
loop calculations. 



3.2.1 Limits on lepton flavor violating interactions 

Using these arguments we can get some indirect bounds on the interactions in eq. @ by using radiative corrections. 
The calculation of the diagram in Fig.[2]remains essentially unchanged if we use as starting point an effective interaction 
with other type of neutrinos or with neutrinos and quarks, thus, if we have an interaction among electrons, u-quarks 

-eP -uP r dP 
-a/3' fc a/3' b af3 

1/3' Lt a/3' Lfc a/3 



or d-quarks with neutrinos with strengths e^, e™^, £^g respectively, we expect a four-fermion interaction among 
electrons, u-quarks or d-quarks with the charged leptons with strengths ce^g, ce™£, ce'S with 



A 



2 In ( — — ) « 0.0027. (46) 



47rs^ \m w 

These interactions give rise to a class of interesting processes like \i~ — > e + e~e~ (BR < 1.0 x 10~ 12 ), r~ — > e + e~e~ 
(BR < 2.9 x 10- 6 ), T~ -> e+e-fi- (BR < 1.7 x 10" 6 ), [iTi -> eTi(T < 4.3 x lO" 12 ), r" -> e-n a (BR < 3.7 x 10" 6 ), 
t~ -> fi-n a (BR < 4.0 x 10~ 6 ), T~ -> e~p°(BR < 2.0 x 10~ 6 ), r~ -> fj,~p°(BR < 6.3 x 10~ 6 ) which are strongly 
bounded from present experiments, thus one can obtain some information on e%^, e^, e^, and similarly for quark 
interactions with neutrinos. We have listed the branching ratios we use to set bounds, so that the limits on the e's 
can be rescaled if the experimental bounds become stronger. 

Consider for example the bound on e*^ from r~ — > [i~ ~e + e~ . The diagram of Fig. 3] with the vertex £g> given by 



eq. © provides an interaction (like that in eq. ([43)1 

C = -2V2G F cef T (eYPe) (Tl 1p Lt) , (47) 



with c given in eq. (|46|) . From this we can immediately compute the branching ratio 

r(r - -> 



BR(t~ -> M^e+e") =BR(t- -> jri/^) Jjr_ _ _ ( = 0.1737 |ce^ | 2 < 1.7 x lO" 6 , (48) 



from where we obtain 

< 0.0031/c< 1.2, 90% CL . (49) 

Note that, as can easily be checked from the Michel parameters, the operators (£^7 p £ a )(e7 p Le) and (Ipj p l a ) 
(ej p Re) give the same contribution to the total decay rate. 

Somehow worse bounds can be obtained for e^, |e^| < 2.9, because the limits on r~ — > e~e + e~ are a bit looser 
and because of the two identical particles in the final state. 

Much more interesting, however, are the bounds that can be set on because of the strong experimental limit 
on BR([i~ — > e~e + e~). In this case the calculation is similar except for a few factors due to the identical particles in 
the final state. We have 

BRQT -> e-e+e") « ^ - e ~ e+e ~) = .5 Ice^f < 1 x 10~ 12 , (50) 
T{fi -> e z^ e ) 



from where we obtain 



e% < 1.4 X 10 _b /c < 5 X 10~ 4 , 90% CL . (51) 



Similar arguments can be used for NSI of quarks with neutrinos. In fact, Fig. [4] also generates interactions like 
eq. (|45p but with the electron fields replaced by u or d-quark fields. These interactions contribute to several hadronic 
decays of the tau lepton. Thus, we can set bounds on the e q J^ (with (3 = e, fi and q = u,d). For instance, assuming no 
unnatural cancellations among the u-quark and the d-quark or the P = L and P = R contributions, we have 



where we included the isospin factor y/2 and used BR(t~ — > v T p ) ~ 0.22. 

From the decay — * e~7r° one obtains similar but slightly worse bounds because the limit on BR(t~ — > e~7r°) 
is worse and because BR(t~ — > v t -k~) is smaller. It is important, however, to remark that decays into p's probe 
the vector channel and decays into 7r's probe the pseudoscalar channel and, in this sense, provide a complementary 
information. 

Using — > yT p° (or r~ — > ^ 7r °) one can a l so se t bounds on e%£. Using decays into p's we obtain \e%E\ < 2.8. 

Finally we can set bounds on e q ^ from [i — e conversion on nuclei. Comparing the experimental upper bound on 
the rate of fi — e conversion to the rate of muon capture, 

R„e = -7^7. — < 4.3 X 10- 12 , 

cr(/U 1 % — ► capture) 

implies (see e.g. [35]) that 

efe & VRi^/c w 7.7 x 10~ 4 . 
We have collected the numerical values of the best bounds in table [3] 



3.2.2 Limits on lepton flavor conserving interactions 

For lepton flavor conserving operators, diagrams like Fig. Q] also give rise to interactions with four charged leptons, 
however the information we can extract from them is not so useful. In addition, in that case, there are other inter- 
esting interactions that can also be generated through radiative corrections. In particular, loop corrections involving 
interactions like the ones in eq. ([2]) can affect in a non-universal way the decay rates of the electroweak gauge bosons 
as shown in Figs. [5][7] This type of interactions also occur in the SM but in the SM they are universal, they have 
the same strength for the different generations. Then, we can bound them by using the various universality tests. 
These are checked, at most, at the level of 0.1%, which is of the order of the SM radiative corrections (apart from 
m t corrections or running a corrections) which are order ^. We expect the one-loop contributions from NSI to give 
corrections therefore we expect, from these processes, to obtain very weak bounds e < 1. Let us see how these 
bounds arise. 
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Figure 5: Contributions of lepton flavor conserving NSI to the vertex of the Z-gauge bosons to electrons, u and 
d- quarks. 

Fig. [5] gives corrections to the hadronic and leptonic decay widths of the Z-boson. The NSI contributions to the 
amplitudes of these processes interfere with the tree-level SM amplitudes and the total correction is proportional to 
the sum of all neutrino couplings. Since e£ P , ef£ are bounded by other means we will consider only v r interactions. 
The calculation of the diagram gives a correction to the couplings of of e,u,d given by (/ = e, u, d) 

lzits w c w \m z j Q=e Air p 
where we have adopted the PDG convention for SM Z-couplings to fermions, 

with gy — gL + gp, gA = 9l — gR (see table [TJ, and (— l) p = +1 for P = L and (— l) p = —1 for P = R. Results for 
g v A are usually presented in terms of an effective g v A at the Z peak which incorporates the SM electroweak radiative 
corrections (but not final state QED corrections). g A is known with a good precision [g A = —0.50111 ± 0.00035 
without assuming universality in the fit) and agrees quite well with the SM prediction for a light Higgs g A w —0.5012, 
therefore additional contributions must be small. Requiring that 6g A is smaller, in absolute value, than the error, and 
taking m(A/m.z) > 1, we find 

J2 £ «« < U7Ts2wC2w 0.00035 x 1.64 a 0.5 , 90%CL . 

P a—e,fi,T 

For the vector part the precision is not so good and the agreement with the standard model is not so perfect: the 
measurement yields gy w —0.03816 ± 0.00047 while the SM prediction for a light Higgs is about gy « —0.037. 
Requiring that the additional corrections are smaller, in absolute value, than the error we obtain 

-|9_„2 2 

Y y (-l) P ef a < w w 0.00047 x 1.64 a 0.7 , 90%CL . 

P a—e,fi,r 

These bounds are comparable to the tree-level limits one can set from e + e~ — * ^vv and, in particular, if we consider 
only couplings to v T and assume no cancellations among different operators we find 

|< P |<0.5. (53) 

If we allow for cancellations among left- and right-handed operators the limit is slightly moved to about 0.6. 

Although the same corrections appear also for light quarks, they only affect the total hadronic decay width of the 
Z , which also contains contributions from s,c and b quarks. One can try to subtract the b and the c contributions but 
the s quark is practically impossible to separate. Therefore, although with additional assumptions (for instance that 
the 6, c and s quark couplings are standard) one can set some bounds on the hadronic e's, the bounds will not be 
comparable to the ones obtained for the electron couplings. 
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Figure 6: Contributions of lepton flavor conserving NSI to the Z invisible decay width. 



On the other hand, Fig. Ogives corrections to the invisible decay width of the Z-boson which, assuming 3 neutrino 
flavours, can be used to set bounds on the e's. In this case non-diagonal neutrino NSI could also contribute. However, 
these contributions do not interfere with tree level amplitudes and are already bounded by lepton flavor violating 
processes. Neglecting the non-diagonal terms, Fig. [6] will be proportional to the sum of all the e's. Evaluation of the 
diagram yields the following correction to the Z-couplings to neutrinos 



\mz J 



U7rs w c w f — ud — 



Assuming that all the other couplings are standard one can obtain information on the above combination of e's 
from the invisible decay width of the Z, Ti nv . SM limits on AiV„, the number of neutrino species, are in fact bounds 
on AT- mv . If only v T NSI contribute to Y\ nv [y e and interactions are already well bounded from other processes) 
we can set bounds on ei£ . Using that the T- mv is proportional to J2 a (9A") 2 ( we assume om y left-handed neutrinos so 
g v £ = <7y° ) one can easily see that 

% = \^N V , (54) 
9a 2 

where in the denominator we have assumed approximate universality in the couplings. Taking only one operator at a 
time and requiring that the NSI contribution is smaller than the error in absolute value (AN V < 0.008) we find (at 
90% CL) 
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Of course for electron couplings to neutrinos we have much better bounds from the leptonic decay widths of the Z. 
For hadronic couplings these are the best limits we have. 

In addition to the processes considered there are other one-loop diagrams that could give interesting contributions. 
For instance, the diagram in Fig. Ogives corrections to the coupling of the W-boson to the electron only. This coupling 
can be absorbed in the definition of Gf affecting equally \x decay and (3 decays (or r — > e v T T> e and 7r — * e D e ). Note, 
however, that it will not affect r — * [iVrV^ or 7r — > \iv^ L . Obviously, these contributions also will affect differently 
W — > eD e and W — * [iv^ or W — > tu t . Therefore, universality limits can be used to set bounds on e%^, the only 
coupling appearing Fig. [7J However, we do not expect interesting bounds since this coupling is already bounded from 
v e e — > v e e scattering. 

Finally one can use the recent determinations of Gf in t decays (such as r — * v T ^,v lll r — * vi^) to set constraints 
[36] on new physics involving the r. Operators of the form (fj p i / T )(f'jpPf) are induced at one loop by the operators 
we are considering (see Fig. [8]). Thus, from the bounds in [36] one can extract limits on the flavour diagonal e{.^. 
However, these limits are loser than the ones already discussed and only affect left-handed couplings. 



Figure 7: Contributions of leptonic flavor conserving NS1 to the vertex of the W gauge boson. 




Figure 8: Contributions of lepton flavor conserving NSI to standard charged current processes. 



In any case, one can use atmospheric neutrino data to set much more stringent limits on these couplings, once we 
use that the are well bounded from scattering experiments. 

4 KamLAND and SNO/SK 

The large mixing angle (LMA) oscillation solution of the solar neutrino problem [37], has been confirmed by the 
anti-neutrino reactor experiment KamLAND. It is a short baseline experiment (matter effects are negligible in the 
present global allowed region by solar and reactor data) measuring electron anti-neutrinos through a charged current 
process. We will assume that KamLAND finds an oscillatory signal, so the neutrino parameters Am 2 oi and Q\i will be 
determined with good precision. Let us stress that the neutral current NSI we are discussing in this work does not affect 
in a significant way the KamLAND observables. On the other hand, the solar neutrino data may be plagued by NSI 
present in the matter potential (evolution in the Sun and in the Earth) and/or by NSI present in the neutrino neutral 
current detection. Therefore the consistency of the KamLAND and the solar neutrino data can give us information 
on NSI that contribute in solar neutrino experiments. In this section, we anticipate the constraints that could be set 
with three years of KamLAND data. 

We will consider solar neutrino data only from SNO and SK (for a discussion on the impact of NSI in the Borexino 
observables, see Ref. [38]). The main source of neutrinos at the energies relevant for SNO and SK is S B neutrinos. 
There is also a small contribution due to hep neutrinos (roughly 0.5% of the total rate if we use the solar standard 
model fluxes) that we consider as a source of systematic error in the experiments. The reason for using only S B 
neutrino data is that the resulting bounds on NSI are almost solar model independent. For simplicity, we will use 
the total number of events measured at these experiments for the daytime and for the nighttime characterized by the 
total rates and the day-night asymmetries. These data are enough to show the bounds that can be reached by this 
set of experiments. A more complete analysis should include spectral SK and SNO information. 

NSI modify the evolution in matter by the effective parameters e' v and e v , described as a function of the funda- 
mental NSI by 03 

e v = J2 (co S e 23 e^-sm0 23 e^)^-, (61) 

f—e,u,d 

e' V = E (4r-(cos 2 e 2 3^ + sin 2 23£ ^-Bin2^))^, (62) 

f—e,u,d 

where Nf is the number density of target particles / in matter, and e v = e L + e R . We are using the fact that the 
3-neutrino evolution in matter can be described in good approximation by an effective 2-neutrino description even for 
e eT ~ 1, because GpN e -C Am 2 tm / ' E v . The neutrino evolution Hamiltonian in the presence of NSI can be found in 

Let us discuss the main effect of NSI in the solar neutrino evolution: 

• e lV modifies the survival probability through the change of the effective matter density seen by the neutrinos. 
If £ is negative, the effective density is smaller and the resonance happens at larger energies, changing the 
recoil energy spectrum and the total rate in SK and SNO. Thus, the evolution of the 8 B neutrinos in the Sun 
is sensitive to negative e' v . On the contrary, a positive e' v increases the effective density and the resonance 
happens at lower energies, having a flat spectrum for the 8 B neutrinos and a slightly affected rate. This would 
imply that there is no bound on positive e' v from 8 J5 neutrinos. However, the large effective density would be 
also in the Earth and it would make an effect in the day-night difference. The measured day-night asymmetry 
(Aun = 2(7V — D)/(N + D) where D(N) are the number of events measured during the daytime (nighttime)) 
bounds positive e' v , 

• e v , if small, modifies the survival probability through the change of the effective mixing 26J. Roughly speaking, 
KamLAND depends on 9 while SNO/SK observables depend on the combinations 6 + e v and the comparison of 
the allowed ranges for the mixing angle extracted from KamLAND data and from the SK-SNO data define the 
e v bounds we can get. 

We consider the following solar neutrino observables: 

II In this analysis we have taken 6*13 = 0, which is a good approximation because solar and KamLAND data are weakly sensitive to #13 
below the CHOOZ bound. 

I H 



SNO Charged Current (CC): NSI appear only in the propagation due to matter effects (only vector couplings 
contribute) , 

[CC] = f B (Pee(e' v ,e v ))cc , (63) 

where 



rhSSM ' 



(64) 



is the 8 B solar neutrino flux normalized to the standard solar model prediction. We denoted [XX] as the observable 
XX normalized to the case of no transition ([XX] = xx(i^ a =i) ) an d ( )xx indicates the observable XX averaged with 
the detector response. P ee is the probability that v e produced in the sun will arrive as v e at the detector. 

SNO Neutral Current: NSI appear mainly in the neutrino detection. At low energies the neutrino-deuteron cross 
section is dominated by the Gamow- Teller transitions, so that the cross section scales as g\, where gA is the coupling 
of the neutrino current to the axial isovector hadronic current E3 [39] 140] 141] 142] [43] . In the SM, gA = g\ — g\ = 1 
and using that the nuclear corrections to gA are the same when the NSI are added we obtain: 



[NC] ~ f B (l + 2e A ), (65) 



where 



e A ~ £ (P ea ) N cKi-e d *i), (66) 

a— e,/^,r 



Thus, this cross section is mainly sensitive to the axial part of the NSI and contains complementary information 
to the oscillation probabilities that depend on the vector part of the NSI. Notice that in the case of absence of axial 
NSI, the NC detection is blind to oscillations (J2 a=e „ T Pea = 1) an d determine the total 8 B flux even if neutrinos 
oscillate |4"4"] . 

SNO-SK Elastic Scattering: NSI appear in the neutrino evolution and also in the neutrino detection (38] . 

[ES] = fB[r e (P ee (e' v ,e v ))ES + 0.157r a (l~(P ee (e' v ,e v )) E s)] , (67) 

where 

re = ^ ^ Ve e{9L + ef e , g R + eg, V |e^| 2 ), (68) 

<JuA9L.gR) frl 



cos 2 0- 



2:i 



0'^e(fl , i)9 , H) 
Vv T e(gL,gR) 



I i ^eL . r-eR, \_eP\2\ 

(Jv T e{gL+e TT ,g R + £ TT ,2_^\e aT \ ). 



(69) 



cr i^ c e(5L + £ ee— ) ^ s gi yen m e Q- ©i the <J Uhl e{gL + e uu— ) an d cr ^e(ffL + £ tt— ) are tne obvious modifications of this 
equation. In the following, we refer to [ES] as the averaged value from SK and SNO electron scattering measurements. 

Regarding the data used in the present analysis, we assume that in the next three years: 

- KamLAND will find a clear signal of oscillations (total rate and energy distortion). From KamLAND, we 
will be able to know the ranges of the solar oscillation parameters, Am^, and 9\2 with good precision. Assuming 
that KamLAND confirm the present best fit point to the solar neutrino data, we expect (P ee (e' v = 0, e v = 0)) = 
0.32(1 ±0.07) [45]. 

- We assume a moderate improvement on the SNO and SK measurements and that the NC is measured indepen- 
dently from the CC in the second and the third SNO phases [46] . The central values correspond to the present central 
measurements : [CC]= 0.35 (1± 0.05), [NC]= 1.01 (1± 0.08) and [ES]= 0.47 (1± 0.03) (the present measurements 
and errors are [CC]= 0.349 (1± 0.057) 0, [NC]= 1.008 (1± 0.125) [2], [ES]= 0.465 (1± 0.032) [47]). 

- We use the day-night asymmetry measured at Super-Kamiokande Adn — 0.021 ± 0.024 [47] . consistent with the 
first SNO measurement A DN = 0.07 ± 0.05 [48] , 

**We are grateful to Jose Bernabeu for an enlightening discussion about this point. 
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- In the calculations, we use 0.4 < cos 2 823 < 0.6, the expected 90% CL allowed range from MINOS [49] . 

Finally, the results presented in this section have been obtained solving numerically the evolution in the Sun and 
in the Earth using the electron and neutron number densities from the standard solar model [50j and the preliminary 
reference earth model [51]. Details on the KamLAND simulation used in this analysis can be found in Ref. [45]. CC 
and NC cross sections were obtained from Ref. [IT]. ES was computed including radiative corrections [52] , 

We present bounds on a given four-fermion operator, assuming that all other four-fermion operators are zero. 
Let us illustrate how the bounds come out for a particular NSI, namely (e"f p Le)(P e "f p Lv e ) . If we consider only this 



operator, eqs. (f6"3"]) - (|U7]) simplify to: 

[CC] = f B {P ee (e£))cc , (70) 

[NC] = J B , (71) 

[ES] = f B (r e (P ee {ef e ))ES + 0.157(1 - {P ee (e^)) E s)) ■ (72) 



In that case, the NC value is a measurement of the 8 B flux, fg= 1.01 (1 ± 0.08). This result for Jb combined with 
eq. ([TO"]) gives {P ee {e%^))cc — 0.33(1 ± 0.09). By inspection, if we compare this averaged rate with the result from the 
KamLAND allowed region, (P ee (e' v = 0,e v = 0)) = 0.32(1 ± 0.07), we can conclude that e%^ must be constrained in 
the range where matter effects from the NSI are small. Using this conclusion in eq. (|T2"|) . e%^ is further constrained 
by the ES dependence on r e . Solving numerically the set of eqs. (|70"]) - (|7"2"]) as a function of fs and and using the 
allowed range of parameters (Am 2 oi ,#i2) from KamLAND, we get the 90%CL bound : 

- 0.08 < ef e < 0.10 . (73) 

More precisely, given the solar observables considered, we computed x 2 so i as a function of the NSI parameter eg, 
the oscillation parameters Am 2 o[ , 612 and #23, and the 8 B flux normalization, fg. We add the x 2 of the KamLAND 
analysis, a function of Am 2 oi and 9x2, to the xLz- Next, we find the marginalized X 2 ( e g) by minimizing the total x 2 , 
Xtotal = Xsoi + XKamLAND> respect to the oscillation parameters and the 8 B flux normalization for each eg. Finally, 
we get the bounds on the NSI parameter by comparing the function Ax 2 (eg) — X 2 ( £ g) — Xmin anc ^ the- statistical x 2 
distribution with 1 dof (Xmin ^ s the minimum of x 2 (total) in the full space of parameters). For the different operators 
(one by one) that appear in the solar observables, we get at 90% CL: 



- 0.2 < e e a L a < 0.3 


(74) 


-0.3 < eg < 0.5 


(75) 


-0.9 < e e a R a < 0.3 


(76) 


-0.25 < e" e L < 0.2 


(77) 


-0.3 < e u a L a < 0.25 


(78) 


-0.2 < elf < 0.25 


(79) 


-0.25 < e£2 < 0.3 


(80) 


-0.2 < eff < 0.25 


(81) 


-0.25 < e d a L a < 0.3 


(82) 


-0.25 < eff < 0.2 


(83) 


-0.3 < 4* < 0.25 


(84) 



where a = h,t and 

- 0.2 < eg < 0.3 (85) 
where / = e,u,d, P = L, R, and a, (3 = e, ji,T with a 7^ (3. 



5 ^Factory 

In this section, we estimate the sensitivity to NSI of ve scattering and i^-DIS at the near detector of a neutrino factory. 
We outline in the next paragraphs the order of magnitude of the limits the near detector could set on e aa (a 7^ t) , 
and e a (a 7^ (3). A more careful analysis follows in subsections 15.11 and 15.21 based on the zz-DIS chapter of the 
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ECFA-CERN Neutrino Factory Study [53] . Mangano et al. [53] discuss various measurements from which sin 2 8w 
could be extracted, with their potential errors. Finally, in section IQ1 we briefly review potential nufactory bounds 
on new charged current neutrino interactions [161 153j . 

For flavour diagonal NSI involving v e and f M , sin 2 Ow measurements at a near detector of a neutrino factory should 
be more sensitive than oscillation probabilities measured at the far detector. e aa will interfere with the SM amplitude 
for v a f — > fct/jO so the diagonal NSI contribute linearly to scattering. They also contribute linearly to the oscillation 
probability at the far detector, via their MSW contribution to the neutrino mass matrix. But there should be more 
events at the near detector, because the near beam is narrow, so a well-instrumented detector, of area larger than the 
beam, can be built at reasonable cost. This suggests that the statistical and systematic errors on e ee and e MM would 
be smaller at a near detector. 

The situation is less clear for flavour changing e a p (a ^ /3) which contribute quadratically to sin 2 9yy . Neglecting 
for the moment systematic errors, one can estimate that the near detector would be sensitive to NSI such that 
^nsi ~ V^Ysm- The number of events in the near detector due to NSI is N^ SI oc e 2 Ng M , where Ng M is the 
number of SM v scattering events. In the analysis of [53] , NJi M ~ 10 8 was taken, so e a p ^ .01 might be seen. At 
the far detector [101 [HI Q~5] , £a/3S ( with a ^ (3) appear linearly in the oscillation probability, because they interfere 
with the SM oscillation amplitude. NSI could disrupt the measurement of oscillation parameters at the far detector 
if nIj SI ^ Af/ sc , where N[ sc is the number of events at the far detector due to oscillations. So e a p ^ sin #13 could 
interfere with the measurement of sin #13 at the far detector |16[ . A recent analysis |54j suggests that a neutrino 
factory can realistically measure sin $13 .01. We find that the near detector would be sensitive to e a p a few x .01, 
so the sensitivities of the near and far detector to flavour-changing NSI are similar. 



5.1 Measuring sin 2 9\y leptonically 



The weak mixing angle could be measured leptonically in the scattering of neutrinos off electrons in the target. The 
error is smaller using the Va, v e from the /i + beam; sin 2 8\y can be determined from 



ov/act = — > ve) + a(v e e — > ve) 



(86) 



with a statistical error of order A sin 2 Ow = 2 x 10 -4 [53 j . To estimate "90 % C.L." limits on NSI, we multiply by 1.6, 
so use A sin 2 6w = 3 x 10 -4 . The second of these two cross-sections, augmented by the contributions of NSI, is given 
in eq. ([4]), and the first is 



crfV^e — * ve) 



2G 2 F m e E v 



\{9l + e e ^f + i£ \ef a \ 2 + (g e R + e^f + £ \ef a \ 



(87) 



Requiring (a u f act - (Tsm) ^ As 2 ^ x dasM/ds w gives 

a a 

+2(1 + gl)et + ]T K L e \ 2 + Z g y£ + i £ K R e\ 

a a 

so the flavour-changing NSI satisfy 



< 6 x 10- 



f;l\ < .04 |ej£| < .02 



\et L A < .02 \e e J}\ < .04 



(89) 

(We neglect e e £ because it is more strongly constrained by fi — > 3e — see section [331 ) The flavour diagonal es satisfy 

|e^|< 0.003 ,\efi\< 0.001 ,\e£\ < 0.0004 , \e**\ < 0.004 (90) 
assuming no cancellations in equation ([88]) . 



5.2 Measuring sin 2 0w in DIS 

In neutrino Deep Inelastic Scattering (DIS), the NC events due to incident v e and v^ cannot be separated. So the 
Paschos-Wolfenstein ratio: 

Rpw = - aNC ^} (91) 

'Tq ^ t here 
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is not available, and s^y cannot be so elegantly disentangled from parton distributions. It was conservatively estimated 
in [53] that a neutrino factory could measure s 2 ^ to one part in 10 3 , via a ratio of the form 



R 



1+r 

where r ~ 0.5 is defined in cq . ([26]) . and the g with lepton subscripts are as defined in equations ([TT]) . (fl"2]) . replacing 
e^/ias required. 

A similar ratio R + can be defined for the beam produced in the decay of a /i + : 



fe) 2 + (9Le) 2 + r[(g L ^) 2 + fee) 2 ] (93) 



1 + r 

R + and R~ can be measured with equal sensitivity, so s 2 ^ can be independently determined from both. The NSI 
appear in one or the other not multiplied by r, so we can set bounds of order A(<7p MjE ) 2 < 9i? ± /9s^ / x 10~ 3 ~ .0005, 
for instance 

tM^BL < 5xl0- 4 (94) 



or equivalently 



e^|<lxl0- 3 , (95) 



and assuming no cancellations among the terms in gpi (eqs.(|lip. (112p ) : 

I^U4e L |<lxl0- 3 , \e%\,\e%\ <9xl0- 4 
\efi\, \e u *\ < 2 x 10- 3 , \e d *\, \e d e ?\ < 5 x lO" 3 

T qL\ \nh\ \ e qR\ \ e qR\ ^ 3 v 1 n~ 2 



e^\e^\,\4^,\47\ < 3x10^ . (96) 



5.3 Charged current interactions 

Non standard neutrino interactions of the form given in equation {2} could contribute at one loop to muon decay or 
charged current v scattering off quarks in the near detector. This would occur via the external dressing of the four 
fcrmion operator with a W loop, as discussed in section 13.21 For instance, exchanging a W between the and e 
legs of e(i> T j p v l j_)(ejpLe) would generate the operator ce(D T j p is e )(ej p Lfi). This would produce a v T from [i decay, 
which could turn into a r in CC scattering off quarks. The sensitivity of the near detector to operators of the form 
2v / 2G_F^(e7p^)(^ / 37 p i' Q ) was estimated in [53] (see also [12l HH [55]) to be 



<r m-4 fx? ~ / / number of A 4 decays 

^ 10 ^ x v^7^v wo (97) 

where N is the number of r events required in the detector for a signal. N ~ 10 was taken in [53]. For c ~ .002, this 
is a weaker bound on e^,e e J^ than from measuring sin 2 9w However, in the opposite limit of c ~ 1 (for instance if 
the NSI are induced by dimension 6 operators), it is clear that flavour changing NSI are more readily detected via CC 
interactions than NC interactions. 



6 Summary 

We have considered non-standard interactions of neutrinos with first generation leptons and quarks, parametrised as 

C eff = -ei P 2V2G F (D alp Lv p )(fYPf) ■ (98) 

We have taken a phenomenological approach, assuming that the new physics which induces these non-standard neutral 
current operators does not generate the SU (2) related charged lepton operators at tree level. This could be the case 
if the operators in ([55]) are of dimension eight or larger. Then the tight bounds on charged lepton NSI do not apply, 
so we have assembled present and future constraints on such operators from purely neutrino processes. 
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We point out though, that even if only neutrino neutral current NSI are present at tree level, they will necessarily 
induce the related vertices with charged leptons at one loop, via W exchange. Moreover, radiative corrections involving 
these neutrino NSI could affect a variety of precision observables. Thus we have also set bounds on the strength of 
the NSI in ([98]) from their one loop effects, which in some cases are more stringent than the tree level ones. 

Our results are summarized in tables [5] and [3] We list the limit on e^Jl that an experiment would set if only one 
NSI operator was present. The limits that arise if cancellations are allowed among the operators are presented in the 
body of the paper. 

Tree level bounds arise from low energy neutrino scattering experiments. We have collected present constraints, 
including the recent NuTeV data, and estimated future limits attainable at the near detector of a neutrino factory 
from measuring sin 2 8yy (both, leptonically and in neutrino DIS). We found that these experiments are more sensitive 
to e{e and ef? than the ^factory far detector. They also provide the best bounds on flavour changing interactions 
involving v T . 

One loop bounds on lepton flavour violating operators are very stringent for . The present experimental limits 
are so strong, that even with the loop suppression, the are constrained to be 0(1O -3 ). The analogous bounds on 
£ ra ( a ~ e > A*) from several tau decays are just order one, although the experimental limits on some of these decays 
could improve, in which case the relevant bound can be rescaled. Regarding flavour diagonal NSI, one loop bounds 
on from the precise measurement of g e A at LEP are of the same order as the tree level limits from e + e~ — > ^vv, 
while the limits on e%^. from the invisible Z width are the only constraints from laboratory data. 

Finally, we have estimated future bounds from comparing KamLAND and solar neutrino data. These could 
constrain at the level of 0(0.3) (which is roughly one order of magnitude improvement of the present limits), and 
set bounds on which are comparable to those from LEP. 

We hoped to show that long-baseline experiments at a neutrino factory would indeed measure oscillation parame- 
ters, since NSI would be seen first in other experiments, like short baseline, high intensity, precision neutrino scattering 
experiments. To what degree have we succeeded? 

Suppose that no evidence for NSI is found; does this mean such interactions can be ignored in long baseline 
oscillation experiments? That is, if the near detector of a neutrino factory sees no NSI, will the far detector measure 
oscillation parameters? The answer we find seems to be that the NSI cannot quite be ignored. The three problematic 
es are e TT , e T)i and e Te . 

The Mangano et al. [53J analysis suggests that the near detector would be sensitive to e Ttl ,e Te <; .02, and the 
Freund et al. analysis [54] says the far detector could see sin6>i 3 & .01. We would like e a p < sin #13 to be sure that the 
es do not confuse the determination of sin #13. Can these two analyses be compared? Freund et al. use more events, 
but the error in the near detector analysis of 53J is systematic, so it is not clear how the errors scale if we increase 
the number of events. The error on the DIS determination of sin 2 Oyy could be significantly smaller than estimated 
in [53], because the source of error is the parton distribution functions, and the near detector will measure these [55]. 
So the near detector might indeed be able to constrain e^,e^ < .01. Also, ^ .01 could possibly be seen at 
ICARUS/OPERA [IT] , where the neutrino beam is largely produced in pion decay, and the detector should be able 
to identify rs. 

e%^,e%^ remain a problem, but could perhaps be disentangled from #13 at a neutrino factory by using the beam 
spectrum, as discussed in |17j . 

Acknowledgments 

We thank Anna Rossi for stimulating our interest in this physics. We also thank Jose Bernabeu, Paolo Gambino, 
Michele Maltoni and Antonio Pich for discussions, and Cecilia Lunardini, Alexander Friedland and Jose Valle for 
comments on the manuscript. This work was partially supported by the Spanish MCyT grants BFM2002-00345, 
BFM2002-00568 and FPA2001-3031, by the TMR network contract HPRN-CT-2000-00148 of the European Union 
and by the Generalitat Valenciana grants CTIDIB/2002/24 and GV01-94. CPG acknowledges support from NSF 
grant No. PHY-0070928. 



on 



Table 2: Flavour conserving four fermion vertices involving two neutrinos and two first generation fermions (ee,dd or uu), 
the best current and the best future 90 % CL limits that can be set on the coefficients 2\/2Gf£ of the four fermion vertices. 
See eq. ([2]) for the definition of e. The limits from processes marked with an asterisk, *>, arise at one loop and are inversely 
proportional to log(A/mw)- We have assumed log(A/myp) > 1 (see section [3~2|l . 



vertex 


current limits 


future limit 




14^1 < 0-5 
(q e A @ LEP)*) 


-0.2 < e% L T < 0.3 
-0.9 < e%* < 0.3 
KamLAND and SNO/SK 




krr 1 < 1-4 

\e^\ < 3 

(r in „)*) 


-0.3 < < 0.25 
-0.25 < < 0.3 
KamLAND and SNO/SK 


(dj p Ld)(i? T jpLi> T ) 


\4r\ < LI 

\ef?\ < 6 


-0.25 < e d T L T < 0.3 
-0.3 < eff < 0.25 
KamLAND and SNO/SK 


(e'y p Pe)(u fl 'y p Lu fl ) 


\ef p \ < 0.03 
CHARM II 


\e%\ < 0.003 
\ef p \ < 0.001 
leptonic at nufact 


(u'y p Pu)(pf 1 j p LUf 1 ) 


\e%\< 0.003 
-0.008 < £«f < 0.003 
NuTeV 


\e%\< 0.001 
\s^\ < 0.002 
in DIS at nufact 


(dyPd)(D plp L^) 


1^1 < 0.003 
-0.008 < e^u < 0.015 
NuTeV 


\e%\ < 0.0009 
|e#f| < 0-005 
in DIS at nufact 




-0.07 < ef e < 0.1 
-1 < e e e f < 0.5 
LSND 


|e^| < 0.0004 
\eg\ < 0.004 
leptonic at nufact 


(uYPu)(D e "f p Lv e ) 


-1 < < 0.3 
-0.4 < elf < 0.7 
CHARM 


|e« e L |< 0.001 
\e^\ < 0.002 
in DIS at nufact 


(dyPd){v elp Lv e ) 


-0.3 < < 0.3 
-0.6 < eif < 0.5 
CHARM 


|el L | < 0.0009 
|eff | < 0.005 
in DIS at nufact 
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Table 3: Flavour changing four fermion vertices involving two neutrinos and two first generation fermions (ee.dd or uu), 
the best current, and the best future 90 % CL limits that can be set on their coefficients 2v / 2Gf£- See eq. §2§ for the 
definition of e. The limits from processes marked with an asterisk, *\ arise at one loop and are inversely proportional to 
log(A/miv)- We have assumed log(A/myp) > 1 (see section [3~2|l . 



vertex 


current limits 


future limit 


(e 7 "Pe)(^ 7p L^) 


K < 1-2 
(r — > /iee)*) 

Kl < o.i 

CHARM II 


K| < 0.04, K| < 0.02 
leptonic at nufact 


{u^ p Pu){D T ^ p Lv^) 


Kl < 2.8 

(t - M* } 
Kl < 0.05 
NuTcV 


Kl < 0.03 
Syy in DIS at nufact 




Kl < 2.8 
(t - 

Kl < 0.05 
NuTcV 


Kl < 0.03 
s^y in DIS at nufact 


(e7 p Pe)(^7 p Li/ e ) 


\ef e \ < 5 x 10- 4 
(P - 3e)*) 




(uyPu)(D^ p Li> e ) 


Kl < 7-7 x 10- 4 
(Ti/x -> Tie)*) 




{dYPd){Dp lp Lv e ) 


K < 7.7 x 10- 4 
(Ti/i -> Tie)*) 






Kl < 2.9 
(r — > eee)*) 
Kl<0.4,l4el <0-7 
LSND 


K| < 0.02, K| < 0.04 
leptonic s^y at nufact 




Kl < 1-6 

(r -> ep)*) 
Kl < 0.5 
CHARM 


Kl < 0.03 
Syy in DIS at nufact 


{dyPd){v Tlp L Ve ) 


Kl < 1-6 
(r -> ep)*) 
K|<0.5 
CHARM 


K| < 0.03 
Syy in DIS at nufact 
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Erratum to Present and future bounds on non-standard neutrino 

interactions 

The renormalisation group mixing between operators that induce Non-Standard neutrino Interactions, such as 
(T7 T 7 p Pii/ M )(e7 P e), and operators that induce flavour changing rare charged lepton decays(such as r — > ^ee), 
is not present at leading log. This is contrary to the claims of Section 3.2 and we thank Mattias Blennow, 
Carla Biggio and Enrique Fernandez-Martinez [T] for drawing our attention to this error. 

The bounds on non-standard neutrino Interactions (NSI) are affected as follows: there are finite contributions 
of the same ordejHl but unlike the leading log, the corrections are renormalisation scheme-dependent. The 
scheme dependence can be cancelled by matching corrections at the high scale, which are model dependent. 
Model and scheme independent bounds could be obtained by performing a NLO calculation. 



Win dimensional regularisation, and unitary gauge. There are additional diagrams in other gauges, but we obtain the same finite part. 



O/l 



7 Details 



Recall that the NSI operators in which we are interested are of dimension 8, are SU(3) x SU(2) x U(l) 
invariant, and are constructed out of four fermion and two Higgs fields. A list of the relevant dimension 8 
operators can be found in [3J. The operators are constructed such that when the Higgs gets a vev, they induce 
the vertex 

(j-ffypalpvp) , (99) 

but not 

(/y/Xw/j) - (frfypatpep) , (ioo) 

where / is one of u, d or en, and a, (3 are flavour indices. To simplify our calculation, we neglect the possibility 
that / = e/,, which implies that W bosons can not be exchanged between the currents. 

We want to obtain the contribution to charged lepton vertices such as Eq. 11001 from the the combination 
of electroweak interactions and the New Physics that induces NSI. We opt to approach this problem from a 
"hard core" Effective Field Theory approach. We imagine the calculation in three steps [3]: 

1. At the mass scale An si of the New Physics inducing NSI, match the full theory onto the Standard 
Model (SM) plus effective operators invariant under SM gauge symmetries. This determines the (model- 
dependent) dimensionless coefficients {Cj} of the effective dimension 8 operators. 

We assume that the coefficients of operators which induce the undesirable vertices (|100p are found to be 
zero to all orders in the new interactions that generate them. This assumption gives the weakest bounds 
on NSI. It needs a significant amount of fine tuning to be arranged; see the discussion in Ref. pQ. It 
could in principle arise due to a symmetry. 

2. Run this effective theory to mw- 

3. Match at mw onto the 5(7(3) x U(l) theory (where there are no Ws); there can be calculable, finite 
terms from the matching at mw- 

The first step we do not consider further. We are interested in the effective theory calculation, which starts 
at step two, with the boundary conditions we have assumed in step one. 

We do not consider quadratic divergences. This is because we work in the Effective Field Theory formalism 
where one calculates in dimensional regularisation, runs with massless RGEs, and at the mass scale of each 
particle, matches the higher scale theory onto the lower scale theory (from which the particle has been 
removed). In this formalism, quadratic divergences appear in the matching, and by assumption in step one 
above, they vanish. 



7.1 Running and matching in broken SU(2) 

To be concrete, suppose that the fundamental theory generates an operator which gives the vertex 
for / = en, with coefficient Cnsi- I n unitary gauge, one can obtain the first vertex of eqn (|100|) by W 
exchange between the neutrino lines (the diagram on the left of figure [9]). It is straightforward to verify, using 
dimensional regularisation, that this diagram is finite. At one loop there is no log, in the lepton massless 
limit, so no mixing during Renormalisation Group running (step two above) of the neutrino vertex of eqn (|99[) 
into the charged lepton vertex of eqn (jlOOp . In our original paper, we incorrectly guessed that this anomalous 
dimension was non-zero so the log was present. However, the diagram has a finite part: at one-loop matching 
at mw, the first vertex Eq. 1 1001 is induced with coefficient 

3m w °nsi hm s 

Unfortunately, this finite part is not scheme-independent. This problem is discussed in Section 1.3. 

We use Naive Dimensional Regularisation for 75. We have checked that we get the same result calculating 
the complete set of diagrams of Fig. [5] in Landau gauge, and in a general R — £ gauge (where the gauge 
dependence vanishes as it should). 



7.2 Running and matching in unbroken SU(2) 

In our EFT formalism, where particles are massless when present in the theory, we could imagine SU(2) to 
be unbroken between An si and mw- This perspective allows a simple, SU(2) invariant reformulation of the 
previous calculation. Consider the operator 

(h p f)Jp SI , (102) 



Figure 9: Relevant diagrams for the broken SU(2) calculation. 



with 

Jp SI = {i a H)~f p {HHp) . (103) 

Recall a, (3 are flavour indices, H is the Higgs doublet, £ a are the lepton doublets and we use the diagonal 
SU(2) contraction 

(HU) = (H* H*_)(^ V \= H*v + H*_e , (104) 
where v and e are lefthanded fields. 

Interactions among four charged leptons will be generated by Eq. 11021 but with Jp 1 replaced by 

J G p hI = (l a H) lp (HH f3 ) , (105) 
with H = ItiH* and iT2 is the antisymmetric SU(2) tensor. 




Figure 10: Current mixing diagrams for an unbroken SU(2). a, b, c, d are SU(2) indices, same for all diagrams. 

We are then interested in the mixing of into the charged lepton interaction current J p LI . We compute 

the one-loop electroweak corrections to this mixing for an unbroken SU(2) theory containing gauge bosons 
and the Higgs doublet. The relevant diagrams are in Fig. 1101 We have used dimensional regularization and 
performed the calculation in a general gauge, checking the cancellation of the dependencies in the gauge 
parameter. The wave function renormalisation does not contribute to the mixing. We find there are no 
logarithmic corrections of NSI into CLI, only a finite term which agrees with the result shown in Sec. 1.1 in 
the theory after SSB. 

Alternatively, we can discuss the mixing between triplet and singlet currents 

J* = HtHl^lf, Jj = h ^ Ah t a r A ^ , (106) 

A 

where {t a } are the SU(2) Pauli matrices. 

NSI and CLI currents are simply related to singlet and triplet currents by 



In this case, after inclusion of wave function renormalisation we find that: 

• Singlet and triplet currents do not mix at all. 

• Singlet and and triplet currents, at one loop, have the same anomalous dimension. 

• In naive dimensional regularisation singlet and triplet currents obtain different finite contributions which 
once translated into the NSI and CLI currents reproduce the finite contribution shown before. 



OK 



7.3 Renormalisation scheme dependence 



If one performs renormalisation group resummation in step two of our EFT prescription, it is well-known 
[3] that one- loop finite corrections are of the same order (in the SU(2) coupling 02) as the two-loop anoma- 
lous dimensions. Therefore, to obtain a renormalisation-scheme independent result, we could compute the 
coefficient of the logarithm arising in the two-loop mixing of J NSI into J s . A useful discussion of the two- 
loop QCD anomalous dimension for a four-quark operator (similar to our four-leg operator, but with scalars 
replaced by fermions), using various prescriptions for 75, can be found in [4]. 

The SU(2) gauge coupling runs more slowly than cxqcd, so one may wonder whether it is neccessary to 
perform RG resummation. Without it, the two-loop anomalous dimension is C(a|) and can be neglected. 
The complete calculation at given order in a-i should be renormalisation scheme independent, if we do 
perturbation theory without resummation (notice that we discuss renormalisation scheme dependence, not 
scale dependence). And, the complete calculation at 0{a.2) includes one loop matching, at the high scale 
An si, of the full theory onto the effective theory. These matching coefficients are clearly model-dependent. 
So if we do not use the EFT prescription of matching-resummation-matching, then the finite bit we have 
obtained must be summed with the (model-dependent) finite bits from the high scale matching. Only if 
we resum, do the scheme-dependent 0(a2(rnw)) terms cancel between the finite part and the anomalous 
dimension. This would allow a result which is formally independent of the high-scale theory. 
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